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The interaction between a superconducting phase qubit and the two-level systems locating inside the Joseph- 
son tunnel barrier is shown to be described by the XY model, which is naturally used to implement the iSWAP 
gate. With this gate, we propose a scheme to efficiently generate genuine multi-qubit entangled states of such 
two-level systems, including multipartite W state and cluster states. In particularly, we show that, with the help 
of the phase qubit, the entanglement witness can be used to efficiently detect the produced genuine multi-qubit 
entangled states. Furthermore, we analyze that the proposed approach for generating multi-qubit entangled 
states can be used in a wide class of candidates for quantum computation. 
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I. INTRODUCTION 

Generation of entangled states of an increasing number 
of qubits has been an important goal and benchmark in the 
field of quantum information''^. Multi-qubit entangled states 
serve as the essential physical resources for measurement- 
based quantum computing^ and quantum error-correcting 
codes^*^. Some of them, such as W stateii, GHZ stated and 
cluster stat9^, have been investigated both theoretically and 
experimentally^iiSiiiiiZiil; however, the experimental prepara- 
tion of multi-particle entanglement has been proved to be ex- 
tremely challenging. To date, the entangled states up to eight 
atoms2 or six photonic qubits''' have been experimentally re- 
ported. As for solid-state systems, due to the difficulty to de- 
couple the qubits with the environments, only the two-qubit 
entanglement of the supercoducting qubits has been demon- 
strated in experiment a'"'-'^-'^ . Therefore, generation of up to 
ten qubits entangled states of the candidates for solid-state 
quantum computation will be a next significant and very chal- 
lenging step towards quantum information processing. 

Among the solid-state systems, superconducting circuit is 
one of the most promising candidates served as hardware im- 
plementation of quantum computer a'^''^-'^ . But the short co- 
herence time limits both of the qubit state manipulation and 
information storage. The loss of quantum coherence in the 
most of solid state qubits is mainly due to the unwanted cou- 
pling of the qubits with the environments. In particular, the 
coherence time would be decay quickly with increasing num- 
ber of qubits because each qubit usually has a control and a 
measurement circuit. Besides fulfilling the manipulations and 
measurements required for the necessary information process- 
ing, all of circuits can also disturb the qubits and lead to deco- 
herence. Therefore, a possible way to experimentally prepare 
more than two-qubit entangled states of solid state systems, 
which have not yet been demonstrated, may need to suppress 
the decoherence from the environments, such as to reduce the 
number of control and measurement lines. 

In this paper, we propose a distinct scheme to prepare gen- 



uine multipartite entangled states for several to ten qubits, 
while the decoherence from the control and measurement 
lines may be minimized by using single control and mea- 
surement setup. The system we have in mind is sev- 
eral to ten of two-level systems (TLSs)^'^i^'-^^-^^i^^-^^-^''i^^--^-'^ 
locating inside a superconducting phase qubit, e. g., a 
current-biased Josephson junction (CBJJ) '^i'^''^ . Recent 
experiment al"'^ i .22,23,24 jj^yg shown that some of TLSs locate 
inside the Josephson tunnel baiTier, while the parameters for 
such TLSs can be detected through spectroscopic measure- 
ments. The lifetime of the TLS is much longer than the deco- 
herence time of the phase qubit, thus TLS can be used as high 
quality quantum memoryai. Furthermore, macroscopic quan- 
tum jump -^"'-^''-^^ has been experimentally demonstrated for a 
hybrid model consisting of a phase qubit and a TLS inside 
the Josephson tunnel barrier^. In particular, it has been pro- 
posed that the TLS itself can be used as qubits for quantum 
computation-^, and typical operations required for the infor- 
mation processing, such as the state initialization, universal 
logical gate operations and readout, have been experimentally 
demonstratec^i. 

Motivated by the progress, in this paper, we show that the 
interaction between the superconducting phase qubit and the 
two-level systems may be described by the XY model, which 
is naturally used to implement the iSWAP gate. With this 
gate, we can effectively generate virous genuine multi-qubit 
entangled states for TLSs, including the W state and mul- 
tipartite cluster states. Moreover, the states of TLSs can be 
manipulated by controlling the interaction between them and 
CBJJ. Such CBJJ-TLS coupling system offers a natural candi- 
date to realize quantum information processing through com- 
bining the advantages of microscopic and macroscopic scale 
systems. In particularly, we found that such entangled states 
may be efficiently detected by measuring the entanglement 
witness^^i^^^ with the help of the phase qubit. Finally, the 
proposed approach to produce genuine multi-qubit entangle- 
ment may be applied to various systems for quantum comput- 
ing, including the trapped ions^ and Cecr^ etc.. 

The paper is organized as follows. In sec. II, we briefly 
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introduce the superconducting phase qubit and TLSs locat- 
ing inside the Josephson tunnel barrier, and then we show 
that such TLSs can serve as qubits for information processing. 
Furthermore, we demonstrate that an iSWAP gate between the 
phase qubit and each of TLSs can be achieved. In Sec. Ill, 
a scheme to implement genuine multi-qubit entanglement of 
such TLSs is proposed. In Sec. IV, the detection of the pre- 
pared multi-qubit entangled states based on the entanglement 
witness is studied. In Sec. V, we show that the proposed ap- 
proach can be used to generate genuine multi-qubit entangled 
states in a wide class of quantum systems, and the paper ends 
with a brief discussion. 



II. TWO-LEVEL-SYSTEMS INSIDE A 
SUPERCONDUCTING PHASE QUBIT 
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The system we consider is a hybrid consisting of a standard 
superconducting phase qubit and several to ten TLSs inside 
the Josephson tunnel barriers, as shown in Fig. 1 . The super- 
conducting phase qubit is a CBJJ, and recent experiments have 
shown that some of TLSs are located inside the Josephson 
tunnel barriers. Furthermore, such TLSs can be considered as 
qubits for the information processing, whereas the Josephson 
phase qubit itself is a 'register' qubit capable of general logic 
operations between TLSs qubits—. The Hamiltonian of the 
phase qubit as shown in Fig. la reads 

where /o is the critical current of the Josephson junction, / is 
the bias current, C is the junction capacitance, $o = ^/2e 
is the flux quantum, Q and 8 are the charge and gauge- 
invariant phase difference across the junction, which obeys 
the convectional quantizing commutation relation [5, Q\ ~ 
lei. For large area junctions, the Josephson coupling ennergy 
Ej = /o$o/27r is much larger than the single charging energy 
Ec = e^/2C. The phase is a well defined macroscopic vari- 
able and quantum behavior can be observed when the bias cur- 
rent is sUghtly smaller than the critical current. In this regime, 
the two lowest energy levels, 1 0) and 1 1 ) , are usually employed 
as two quantum states to form a so-called phase qubit. Trun- 
cating the full Hilbert space of the junction to the qubit sub- 
space, the phase qubit Hamiltonian can be written as. 

Hp = -iwiocr^. (1) 

where wio is the frequency difference between |0) and 
The states of the qubit can be fully controlled with the bias 
current in the form of 

lit) = Idc + Iif{t) + Ifiwcit) cosujiQt + If_,ws{i) sinwioi, 
where the classical bias current is parameterized by I^wc, 

The transition frequency loqi can be measured using 
spectroscopy^"* and is a continuous function of bias current. 
It was found that some TLSs may locate inside the Josephson 



FIG. 1: (Color online) The phase qubit and TLSs locating inside the 
Josephson junction tunnel barrier, (a) Schematic show of the hy- 
brid system. The system is a CBJJ coupled with some embedded 
TLSs. (b) Schematic energy level diagram for a junction coupled to 
a TLS. The ground state (the excited state) of the CBJJ is denoted as 
|0) (|1)), and the frequency difference is ujuy, \g) and |e) represent 
the ground state and the excited state of TLS with level spacing ujr. 
(c)Schematic show that Qubit frequency Lj\o/2-n vs current bias for 
the capacitively shunted design tunnel junction. Similar to the exper- 
imental observation, we plotted about ten splittings observed in the 
spectroscopy. 



tunnel barrier (see Fig. 1(a)). A TLS is understood to be an 
atom, or a small group of atoms, that tunnels between two lat- 
tice configurations^i^. A TLS can induce an energy splitting 
in the curve of the function looi{I) , as shown in Fig. 1(c). 
The splittings Aj e [20, 100] MHz were observed and they 
are separated by 5f ~ 200 MHz on average-- in the spec- 
troscopy. Therefore, one can characterize the positions and 
sizes of TLS in the CBJJ energy spectrum through spectro- 
scopic measurements. When the register qubit is detuned from 
the TLS by 6f, the effective coupling strength is described by 
A| /4(5/. Therefore, with a splitting magnitude of several tens 
of MHz, it is reasonable to assume that only a single TLS sat- 
isfies the near-resonance condition while the other TLSs are 
far off -resonance. Furthermore, as the number of TLS mainly 
depend on the tunnel junction area, the design with small area 
and external low-loss capacitor (keeping the critical current 
constant) will greatly improve its performance^. It is now 
possible to obtain a tunnel junction within about ten useful 
TLSs in the CBJJ spectroscopy ranging over ~ 2 GHz with 
the improved design2&. 

As shown in Fig. 1(a), each TLS can be modeled as a 
charged particle that can tunnel between two nearby different 
positions with different wave functions \R) and \L) (corre- 
spond to critical currents and I^) within the tunnel barrier. 
The interaction Hamiltonian between the resonators and the 
critical-current is^S 

H^nt = -^^cosd<E)\R){R\-^^cos6<»\L){L\. (2) 
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Assume a symmetric potential with energy separated by huil 
for the ith TLS, then the ground and excited states are \g) ~ 
{\R) + |L))/\/2 and |e) ~ {\R) - \L))/V2. In practical 
experiments, the junction is biased near its critical current, 
(5 ^ (5 + 7r/2 with (5 <C 1, thus cos (5 S. In the basis 
{\0g), \lg), |0e), |le)}, as shown in Fig. 1(b), denoting Sij = 

{i\S\j), usually Si-, 



and Sn 



2tt 



Then 



the Hamiltonian of CBJJ-TLS system becomes^^ 



total 



E 

i 



(3) 



where the effective coupling strength between the CBJJ and 

thej-thTLS is 5, = ^Jl^^f^ 

al ^ operate on the jth TLS states. The CBJJ and the jth 
TLS are tuned into resonance when the coupling term satisfy 
jftwio — huj^l < Sj. The parameters uj^ and Sj of the jth 
TLS are usually unknown, but they can be independently de- 
termined by the energy splitting22i2^ on the spectroscopy, as 
shown in Fig. 1 . 

In the interaction picture and under the rotating-wave ap- 
proximation, the effective interaction between the jth TLS 
and the phase qubit is given by 



tary operator described in Eq.® by controlUng the interaction 
time between CBJJ and TLSs. 

We assume that there are N TLSs locating inside the 
Josephson tunnel barrier, while the coupling constant Sj ( j = 
1 , 2 • • • , N) have been measured by using the spectroscopy. 
We first show that the N-qubit W state can be achieved simply 
by two steps. 

The first step is to initialize each TLS qubit in any general 
state \tp)j = aj\g)j + Pj\e)j to the ground state \g)j by per- 
forming an iSWAP operation between this TLS and the phase 
qubit. We have demonstrated in the previous section that such 
iSWAP gate can be realized by switching on the interaction 
between the jth TLS and the phase qubit with the fixed time 
Tj = Tr/2Sj. The ground state of the jth TLS is initialized by 



The Pauli matrices the following transformation 



(4) 



So the effective CBJJ-TLS interaction is described as the XY 
model. Assume that the interaction between CBJJ and the 
jth TLS is tuned to the resonance with time t, one obtain the 
evolution operator given by 



Uj {t) = exp 



(5) 



between them. This operator results in an oscillation between 
\lg) and |0e) at a frequency Sj, i.e.. 



\0g)^\0g),\le)^\le), 

\lg) cos{Sjt)\lg) ~ ism{Sjt)\Oe) , 

|0e) cos{Sjt)\Oe) - ism{Sjt)\lg). 



(6) 



Note that the iSWAP gate between CBJJ and TLS can be ob- 
tained when tj = tt/ {2Sj) = Tj. In the paper, the TLS states, 
\g) and |e), represent our logic qubit. The experiments have 
shown that the lifetime of such qubits is sufficiently long to 
carry out precise gate operations between the phase qubit and 
TLSs. In addition, we assume that the states of different TLSs 
are well separated from each other in frequency. Therefore, 
by adjusting the bias currents, the phase qubit and a TLS can 
be tuned into and out of resonance (turning on and off their 
coupling), i.e., it allows the independent manipulation of each 
TLS. 



in. IMPLEMENTING THE W AND CLUSTER STATES OF 
TWO-LEVLE-SYSTEMS 

We now turn to demonstrate that the genuine entanglement 
of several to ten TLSs could be generated through the uni- 



[tSWAP]^^^P^\^P), ® |0) \g), <E> (a,|0) - z/3,|l» (7) 

with [iSWAP](j p'f denoting the iSWAP operation between 
the jth TLS and the phase qubit, provided that the initial state 
of the phase qubit is |0). The ground state <S>j'=i \9)j for 
TLSs is then realized after an iSWAP gate is performed be- 
tween each TLS and the phase qubit, hybrid an operation to 
initialize the phase qubit to the state |0) between two iSWAP 
gates. 

In the second step, the phase qubit is first tuned far off- 
resonance with all TLSs and excited on the |1) state. Then 
the phase qubit is adiabatically tuned into resonance with one 
of TLSs (f.g., the jth TLS), effectively turning on the cou- 
pling between this TLS and the phase qubit with the time 
tj. Sequentially, the register qubit is interacted with every 
TLSs only once with appropriate time tj, and then an oper- 
ator Uj = e~*^int*j/'i between the phase qubit and the jth 
TLS is performed. In this case (as shown in the Fig. 2(a)), the 
final state of the TLSs and the phase qubit becomes 

N N 

1=1 j=i 

N N 

n C^; (cos(5iii) 1 1) 1.9) 1 - * sin(5iii) |0) |e) i) (g) 

1=2 j=2 
N 

^[]C/acos(5iti)|.g)i(cos(^2t2)|l)|.g)2 



1=3 



N 



-isin(52t2)|0)|e)2)-isin(5iii)|e)i|g)2|0)](g)|5), 

j=3 



N /-I N 

i\0)Y,I[^°<Sjt,)sm{Siti)(^ 

1 = 1 j = l k^l 



\9)k\e)i. (8) 



If the interaction time between the jth TLS and the phase 
qubit is chosen specifically as tj = arcsin(^^^^^^=j), i.e., 

Hj^i cos{Sjtj) sin(S';t;) = then the final state of Eq.© 
indeed becomes the standard W-state of the N TLSs described 
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By using the above process, two important entangled states 
may be realized. Firstly, the Bell state of two arbitrary TLSs 
j and k can be achieved as 

U,Uk\l)\9),\9)k^ 7|(l3>.|e)fc + |e),|.9)fc)|0>, (9) 

when the interaction time are chosen as tj = Tj/2 and tk — 
Tfe. Secondly, the W state of arbitrary three TLSs, such as j, k 
and l-th TLS, may also be obtained through 

U,U,Ui\l)\g),\g),\g)i^ 



(a) 



x/3 



{\9)j\9)k\e)i + \9),\e)k\9)i + \e) ,\g),\g)i)mm 



Tj/3, tk = Tk/2, and 



state 
N TLSs 



\Cn) 
may be 



by choosing the interaction time tj 
ti = T(/2, respectively. 

Furthermore, the cluster 

plemented by a similar process. The process is also two 
steps, (i) The first step is to initialize the A^th TLS to 
the ground state \g)N and all the other TLSs to the states 
= -j=[\g)j + \e)j), while initialize the phase qubit to 

1+) p = -i=(|0) + |1)). (ii) Sequentially perform the iSWAP 
operations between CBJJ and each of TLSs. Note that, to 
cancel the single-qubit phase factor, a Zp ( which denotes a 
tt/2 rotation around the z axis) pulse needs to be applied to 
the phase qubit both before and after each iSWAP operation. 
After that, all N TLSs may be connected to make a large 
cluster chain (as shown in the Fig. 2(b)), i.e., 

[iswap](^N.p) n {^pi^^^M{j.P)[Z]p} 
i=i 

{(^\+)j}®\9)n\0) 

N 



1 



(|.g),a^+i + |e),)|0) 



\Cn)\0). 



(11) 



Moreover, The chain cluster states can be connected to pro- 
duce higher dimensions cluster states by repeating iSWAPs^. 



IV. DETECTION OF THE GENUINE MULTI-QUBIT 

ENETANGLEMENT 



To ensure the obtained state is the desired multi-qubit en- 
tangled state, one must detect the state of TLSs. Because 
the TLSs cannot be directly measured, one must previously 
perform iSWAP operation to transform the state of the jth 
TLS to the phase qubit (up to a correctable Z rotation), then 
the information can be read out through measuring the phase 
qubit with quantum state tomography (QST)?-. Currently the 
measurement?' of the phase qubit with high-fidelity {F = 
0.96) can be completed with a short time (less than 5 ns). The 
readout technique is achieved by applying a short bias current 
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FIG. 2: Quantum circuits for generating (a) W state and (b) cluster 
state of N-TLSs. 



pulse dl{t) that adiabatically reduces the well depth AU /fiiOp, 
so that the first excited state lies very near the top of the well 
when the current pulse is at its maximum value. In this way, 
one can read out the states of TLSs one by one. 

For detecting two-qubit entangled states, one can use tomo- 
graphic state analysi s'"''"^-'^ to reconstruct the density matrix 
with near-unity detection efficiency. This is achieved by single 
qubit rotations and subsequent projective measurements. For 
the two-qubit system, a convenient set of operators is given by 



16 operators ci^'' ( 



("21 

fjg (a, /3 = 0, X, y, z), where denote 



Pauli matrices of qubit j. The reconstruction of the density 
matrix p is accomplished by measuring the expectation values 



.(1) 



^0 



It has been shown that only the expectation 
of az in the phase qubit can be measured; however, the other 
direction measurements can be achieved by applying a trans- 
formation that maps the detected eigenvector onto the eigen- 
vector of (Jz before detecting. To obtain all 16 expectation 
values, nine different settings have to be used. 

A disadvantage of the tomography is that the operators re- 
quired to detect the entanglement are growing exponentially 
with the number of qubits. However, if one knows about 
some priori information about the generated entangled state, 
one can use entanglement witness operator to distinguish and 
characterize the -partite entangled stata^^i^i^. The W state 
and cluster states proposed in the paper are actually two typ- 
ical types of genuine multipartite entangled states (i.e., the 
state whose reduced density operator of any subsystem has 
rank larger than 1). The entanglement witness W is an opera- 
tor such that for every product state 

Tr(pW) > (p e Ss) 

with Ss denoting the set of separable states. From the defini- 
tion of the operator W, it is clear that the witness has a posi- 
tive or zero expectation value for all separable states, and thus 
a negative expectation value signals the presence of genuine 
multipartite entanglement. In order to measure the witness W 
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of the generating entanglement of the TLSs proposed here, we 
should decompose the witness operator into a sum of locally 
measurable operators. By appropriately constructing W, the 
required measurement settings are much less than the require- 
ment of the tomography. 

As for the 7V-qubit W state denoted as | Wn ) , one of the uni- 
versal methods to construct the entanglement witness Wwn is 
given by 



Wn 



N -1 
N 



\Wn){Wn\ 



with / denoting the identity operator Especially, it has been 
proven that the optimal decomposition of the witness Wws is 
given by^ 



m){W3\ 



— [17- 
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I<7z<Jz) 



I<7zI + II(7z) 



«i3 _ j-j 



^z 
Oz 



(12) 



This decomposition requires five measurement settings, 
namely af^ and ((cr^ + (t^)/-\/2)®^, rj = x,y. Since three 
qubits entangled states have not been realized in solid state 
systems, the above optimal decomposition is very useful to 
detect a purely genuine entangled state in near future exper- 
iments. In addition, a universal method to construct the wit- 
ness Wwn which requires N"^ — iV + 1 measurement settings 
is developed in Refj^. 

We now turn to address a very efficient method to construct 
the entanglement witness Wcn of '^he 7V-qubit cluster state 
\Cn)- By using the stabilizing operators S) " of the cluster 



state (i.e., S", 



the required settings for the entanglement witness increase 
at most polynomial with the number of qubits, and specially 
for the cluster state, only two local measurement settings are 
needed for any number of qubits. 



V. GENERALIZATION AND CONCLUSION 

Actually the proposed approach to produce genuine multi- 
qubit entanglement may be applied to a wide class of the 
candidates for quantum computation, examples including 
trapped-ion quantum computation and quantum computation 
based on Cgo etc.. For concreteness, we generalize this 
method to achieve multi-qubit entanglement of trapped ions 
proposed in Ref.^^. Trapped atomic ions remain one of the 
most attractive candidates for the realization of a quantum 
computer, owing to their long-lived internal qubit coherence. 
The central challenge now is to scale up the number of trapped 
ion qubits. However, scaling the ion trap to interesting num- 
bers of ions poses significant difficulties. In Refj^, the au- 
thors proposed an interesting scheme to scale the ion qubits. 
As shown in Fig. 3, ion qubits locate in a two-dimensional ar- 
ray of micro-traps, where the distance between the ions in the 
plane can be very large, since no direct interaction between 
them is required. A different ion named the head ion can 
move above the plane of the ion array. By switching on a laser 
propagating in the perpendicular direction to the plane, one 
can perform the two-qubit gate between the target ion and the 
head ion. In particular, if the two-qubit iSWAP gate described 
in the present paper can be performed between the head ion 
and all ion qubits in the array, one can produce a large number 
of the cluster state or W state. IF the cluster state of the ion 
qubits in this two-dimensional array can be achieved, an one- 
way quantum computation may be implemented, since single 
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(Cn) 



" X " Z 1 



s 



(13) 

rp-i)ai^V(^+i) (j = 2,3,...,7V-l),(14) 



,(Af-l)^(JV) 



(15) 



one can construct the entanglement witness Wcn ' which de- 
tects genuine A^-qubit entanglement around the A^-qubit clus- 
ter state, given by^ 



Wc„ =37-2 



n 



s. 



{C^ 



n 

odd k 



(Cn) 



(16) 



A remarkable feature of this entanglement witness is that only 
two local measurement settings, i.e.. 



motion 




target 



^(2) . , 
(2) . , 



^ X 



all) ® al^) . . . al^-i) ® a^^ ■ ■ ■ al^-D ® ai^' 



X 



UN) 



are needed independent of the number of qubits. Comparing 
with quantum state tomography that the number of measur- 
ing settings increase exponentially with the number of qubits. 



FIG. 3: Scalable quantum computer: two-dimensional amy of 
micro-traps. 
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qubit gate for each ion qubits and the measurement operators 
can be easily reaUzed for trapped ions. 

Furthermore, if each trapped ions in Fig. 3 is replaced by 
a qubit consisting of Cgo'^^, the genuine multi-qubit entan- 
glement for such qubits can also been achieved by using the 
scheme we proposed here. Noted that it is difficult to realize 
the strong couplings between all nearest neighbor qubits. So 
the approach proposed here is promising. 

In conclusion, we have presented an efficient scheme to 
realize and detect the multi-qubit entangled states of TLSs 
locating inside the Josephson phase qubit, and the proposed 



method can be able to applied in a wide class of candidates 
for quantum computation. 
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